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ABSTRACT 

The pulse shapes detected during multiple outbursts of SAX J1808. 4-3658 are analyzed in order to 
constrain the neutron star's mass and radius. We use a hot-spot model with a small scattered-light 
component to jointly fit data from two different epochs, under the restriction that the star's mass 
and radius and the binary's inclination do not change from epoch to epoch. All other parameters 
describing the spot location, emissivity, and relative fractions of blackbody to Comptonized radiation 
are allowed to vary with time. The joint fit of data from the 1998 "slow decay" and the 2002 "end 
of outburst maximum" epochs using the constraint i < 90° leads to the 3er confidence constraint on 
the neutron star mass 0.8M Q < M < 1.7M and equatorial radius 5km < R < 13km. Inclinations as 
low as 41° are allowed. The best- fit models with M > 1.0M Q from joint fits of the 1998 data with 
data from other epochs of the 2002 and 2005 outbursts also fall within the same 3ct confidence region. 
This 3(7 confide nce region allows a wide variety of hadronic equations of state, in contrast with an 
earlier analysis (jLeahv et al.l 120081 ) of only the 1998 outburst data that only allowed for extremely 
small stars. 

Subject headings: equation of state — pulsars: individual: SAX J1808. 4-3658 — gravitation — stars: 
neutron — stars: rotation — X-rays: binaries 



1. INTRODUCTION 

The discovery ()Wiinands fc van der Klislll998f ) of the 
first accretion powered millisecond pulsar, SAX J1808.4 
3658 (hereafter SAX J1808)), confirms the proposal that 
neutron stars in low-mass X-ray binaries are the progen- 
itors of millisecond-period pulsars. SAX J1808 is now 
one of eleven known accre ting millis e cond X-ray pulsars 
(see lWiinandi (|2006l) and iPoutaner] (|2006t ) for reviews . 
The X-ray pulsations are most likely produced from ac- 
cretion onto the neutron star' s mag netic poles (see, e.g., 
Fig. 12 of IGierlihski et al.l (|2002f lL Spectral models 
(jGierlinski et al.1 120021 ) provide strong evidence that the 
X-rays correspond to blackbody emission from a spot on 
the star which is then Compton scattered by electrons 
above the hot spot. Since the pulsed light is emitted from 
(or very close to) the neutron star's surface, accreting 
millisecond X-ray pulsars are excellent targets for light- 
curve fitting in order to constrain the neutron star equa- 
tion of state (EOS). The X-ray light curve depends on 
the intrinsic properties of the emission spot (size, shape, 
location, and emissivity), as well as the neutron star's 
properties (mass, radius, and spin). Constraints on the 
star's mass and radius lead to constraints on the EOS of 
supernuclear density material. 

The first pulse shape analysis (jPoutanen fc Gierlihskil 
(l200l . hereafter PG03) for SAX J1808 provided inter- 
esting constraints on the neutron star's mass and radius. 
However, this analysis did not take into account two 
effects that are important for rapidly rotating neutron 
stars: var iable time delays due to light travel time across 
the star ("Ca deau et al. 20051) and the oblate shape of 
the star ()Cadeau et al.l 20071 ). The pulse shape analysis 
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including these effects was done bv iLeahv et al.l ((2008), 
resulting in smaller mass and radius than found in PG03. 
Since the pulse shapes of SAX J 1808 were variable during 
the time period analyzed by PG03 and during subsequent 
outb ursts by SAX J 1808 (s ee the comprehensive analy- 
sis bv Hartma n et all ((2008)), it is important include the 
time variability in the analysis. In effect, the results of 
an analysis of a time-averaged pulse profile and an anal- 
ysis of individually different pulse profiles are expected 
to be different. Here we carry out an analysis of a set of 
different pulse profiles from different time periods from 
SAX J1808. In addition, fitting multiple pulse profiles 
with the same model should provide a much more strin- 
gent test of the validity of the model than fitting of a 
single pulse profile. 

The outline of our paper is as follows. In section 2, 
we introduce the different pulse shapes observed from 
SAX J1808, which come from different phases of the 
1998, 2002 and 2005 outbursts. In section 3, the hot 
spot model is described and we note that an additional 
scattered-light component is required in order to fit all 
of the pulse shapes. In section 4, the results are given 
for joint fits of the pulse shape data, including resulting 
mass, radius and inclination values. We conclude with a 
general discussion of the results, comparing with previ- 
ous work. 

2. DATA 

SAX J1808 was observed by RXTE during the 1998, 
2002, 2005 and 2008 outbursts. In this paper we concen- 
trate on the publicly available data from 1998, 2002 and 
2005. Since the pulse shapes detected during the 2008 
outburst are very sim ilar to the pulse shap es detected in 
the earlier outbursts (jHartman et al.l[2009f ) , the addition 
of the 2008 data is not e xpected to add signifi cant in- 
formation to the analysis. lHartman et al.l (|2008l) showed 
that the pulse shape is variable during an outburst and 
found shorter time periods over which the pulse shape is 
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stable, as sho wn in Figure 3 o f their paper. The notation 
introduced bv lHartman et al.l (|2008l) for each time period 
of an outburst is: 1, burst rise; 2, beginning of outburst 
maximum; 3, end of the maximum; 4, slow decay; 5, 
steep luminosity drop; and 6, the flaring tail. We make 
us e of the same t i me in tervals and notation presented 
bv lHartman et al.l (120081 ) . For instance, 1998B4 refers to 
the 1998 outbur st during slow decay st age, shown in Box 
4 of Figure 3 in lHartman et al.l (|2008l ). In our analysis, 
we can only make use of data that has a high signal- 
to-noise ratio so that meaningful fits to the data can be 
performed. This results in a total of 7 time intervals, 
which are used in this paper: 1998B4, 2002B3, 2002B4, 
2005B1, 2005B2, 2005B3, and 2005B4. In the case of 
2005B4 the data from 2005B4A and 2005B4B were com- 
bined in order to provide a stronger signal. The flux from 
the source is very low during stages 5 and 6 so unfortu- 
nately useful pulse shapes are not available. In all cases, 
the pulse shapes are binned to 32 bins per pulse period. 

The data for each time interval are binned into two 
narrow energy bands, a low-energy band from 3-4 keV 
and a high-energy band from 9-20 keV. This choice 
of e nergy bands is m otivat ed by previous spectral mod- 
els (|Gierlihski et al.ll200l IPoutanen fc Gierlinskil l200l 
for the pulsed emission that includes two components: 
blackbody emission and a Comptonized powerlaw. The 
high-energy band is dominated by the powerlaw compo- 
nent, while the blackbody component makes a signifi- 
cant contribution to the low-energy band. This choice 
of bands avoids the energy region between 5 and 8 keV 
which includes the iron line contributio n from the disk 
(jCackett et al.ll200l IPapitto et al.ll2009l ). 

In this paper we perform joint fits between different 
time periods. In order to provide the most useful con- 
straints on the model, we need to identify the time inter- 
vals in which the pulse shapes are least like each other. 
To this end, we compute reduced % 2 values comparing 
each of the 21 possible pairs of time intervals using both 
energy bands described above with relative normaliza- 
tions a free parameter. The lowest reduced \ 2 i s 1-17 
(62 degrees of freedom) for the pair 2002B3 and 2005B3 
showing that the pulse shapes are not statistically differ- 
ent from each other. The next lowest reduced x is 3-96 
which has chance probability < 10~ 16 . Thus all other 
pairs of pulse shapes are significantly different. For this 
reason, we omit the 2005B3 pulse shape from our anal- 
ysis since it has larger error bars that the 2002B3 pulse 
shape. The data for the two energy bands for each of 
these time periods are shown in Figures Q] - [SJ The low- 
energy band data are shown using squares in Figures [T] 
- [6] while circles denote the high-energy band data. The 
largest reduced \ 2 values result from comparisons of the 
1998B4 data with any of the other data, ranging from 
a low of 22.5 (1998B4 and 2005B1) to a high of 63.0 
(1998B4 and 2005B4). The large differences can be eas- 
ily seen by eye: the pulse modulation is about 10% for 
the 1998B4 data while it is only about 5% for all of the 
other time periods. Thus we find that the 1998 outburst 
pulse shape is significantly different from all of the pulse 
shapes resulting from the 2002 and 2005 outbursts. 
We now compare the data used in this pa- 
with the data used in previous ana lyses 
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Figure 1. Data from the 1998B4 period of the 1998 outburst 
(squares 3-4 keV, circles 9-20 keV). The solid and dashed curves 
correspond to the best-fit model that results from a joint fit of the 
1998B4 and 2002B3 data sets. This model (see row 1 of Table 
has 2M/R = 0.6, M = 1.31M© and \ 2 = 116 for 110 degrees of 
freedom. 




Figure 2. Data from the 2002B3 period of the 2002 outburst 
(squares 3-4 keV, circles 9-20 keV). The solid and dashed curves 
correspond to the best-fit model that results from a joint fit of the 
1998B4 and 2002B3 data sets. This model (see row 1 of Table EJ 
has 2M/R = 0.6, M = 1.31M Q and x 2 = H6 for 110 degrees of 
freedom. 
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Figure 3. Data from the 2002B4 period of the 2002 outburst 
(squares 3-4 keV, circles 9-20 keV). The solid and dashed curves 
correspond to the best-fit model that results from a joint fit of the 
1998B4 and 2002B4 data sets. This model (see row 1 of Table 
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of freedom. The best-fit light curves for the 1998B4 data are not 
plotted since the curves are not easily distinguishable from the 
curves plotted in Figure [T] 
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Figure 4. Data from the 2005B1 period of the 2005 outburst 
(squares 3-4 keV, circles 9-20 keV). The solid and dashed curves 
correspond to the best-fit model that results from a joint fit of the 
1998B4 and 2005B1 data sets. This model (see row 1 of Table [6j 
has 2M/R = 0.6, M = 1.36Af and \ 2 = HI HO degrees of 
freedom. 




Figure 5. Data from the 2005B2 period of the 2005 outburst 
(squares 3-4 keV, circles 9-20 keV). The solid and dashed curves 
correspond to the best-fit model that results from a joint fit of the 
1998B4 and 2005B2 data sets. This model (see row 1 of Table 
has 2M/R = 0.6, M = 1.2OAf and \ 2 = 109 for 110 degrees of 
freedom. 
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Figure 6. Data from the 2005B4 period of the 2005 outburst 
(squares 3-4 keV, circles 9-20 keV). The solid and dashed curves 
correspond to the best-fit model that results from a joint fit of the 
1998B4 and 2005B4 data sets. This model (see row 1 of Table EJ 
has 2M/R = 0.6, M = 1.19M and x 2 = 109 for 110 degrees of 
freedom. 



from the 1998 outburst that approximately spans 
the 1998B4 and 1998B5 time periods, binned into 
16 time bins per period, and separated into a 3-4 
keV low-energy band a nd a 12-18 keV h igh-energy 
band. I n the analysis bv iLeahy et all (|2008f ) the same 
data as iPoutanen fc Gierlihskil (|2003l ) was used, but 
a bolo metric light curve constructed by iPapitto et al.l 
(120051) using data from a very short period (22 hours) 
of the 1998B4 data was also used. The constraints on 
the mass and radius arising from the bolometric light 
curve were not very strong compared to th e constraints 
arising from the two-band light curves (see ILeahy et al.l 
(2008)). In the present analysis, we require time periods 
over which the pulse shape is s t able, thus we use those 
determined by Irlartman et al.l (|2008l) . We use pulse 
shapes with 32 time-bins per period because they show 
some features that are smoothed out when binned down 
to 16 time-bins per period. The high-energy band that 
we use in the present analysis is somewhat wider than 
the high-energy band used by ILeahy et al.l ((2008) . We 
have confirmed that the change in the high-energy band 
does not significantly change the best-fit models of the 
1998 outburst data. For the pulse shape analysis we 
added systematic errors equal to the statistical errors, 
as is common practice. 

3. HOT SPOT MODEL 

The models used in this analysis are similar to those 
used bv ILeahy et all (120081 ) to model SAX J1808 and by 
ILeahy et all (j2009D to model X TE J1814-338. (These 
model s are similar to that of IPoutanen fc Gierlihskil 
(2003) for SAX J1808, but include ti me-delays and 
oblate ness, which are important, e.g. see iMorsink et all 
(|2007t )). The model for each pulse shape is composed of 
three components: (1) Comptonized flux in the high- 
energy band, (2) Comptonized flux in the low-energy 
band, and (3) Blackbody flux in the low energy band. 
For each component, the observed flux, integrated over 
the energy band, is given by 



^ =W +i '(l-a lM ) 



(1) 



where the subscript "i" takes on values of 1, 2, or 3 for 
the 3 components, Ii is a normalization, 77 is the Doppler 
boost factor that depends on the pulse phase, Ti is the 
effective power- law index for the component, a.; is the 
anisotropy parameter, and fj, is the cosine of the angle 
between the direction that the photon is emitted and 
the normal to the surface at the location of the emitter. 
As in previous models, we take 0,3 = for the blackbody 
component, and set a± — a 2 = a for the Comptonized 
components. The spectral model fixes the values of the 
power-law indices 1^ as well as the relative normalization 
of the low-energy band components through the param- 
eter b — F 3 /F 2 . In these models, it is assumed that 
the Comptonized and Blackbody emissions both origi- 
nate from the same region on the surface of the star. 

The flux in each time-bin is com puted through 
the o blate Schwarzschild approximation (Morsink et al. 
120071 ) which accounts for the light-bending caused by 
gravity, the time-delays caused by the finite travel time 
for photons across the star and the oblate shape of the ro- 
tating star. Since the light-bending depends only on the 
ratio of the neutron star's mass to radius at the spot lo- 
cation, the star's gravitational radius, 2GM/Rc 2 , is kept 
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fixed for a series of fits. (For simplicity, when referring to 
the gravitational radius, we use units with G = c = 1.) 
The lowest x 2 f° r a given 2M/R is found. This allows 
the most efficient use of "look-up tables" for the light- 
bending. Once the the ratio of 2M/R has been specified, 
there are eight free parameters for a fit to one pair of 
2-energy band pulse shapes (say 1998B4): R, 9, i, I\, I2, 
b, a, and <f>. The angle 4> is an arbitrary phase. The angle 
between the spot's centre and the spin axis is 9, and the 
observer's inclination angle is i. We assume that the or- 
bital and spin angular momenta are aligned, so that the 
inclination angle i is the same as the binary's inclination 
angle. Since the angles 9 and i are both measured from 
the "north" spin axis, an angle smaller than 90° corre- 
sponds to a spot or observer in the northern hemisphere, 
while an angle larger than 90° corresponds to the south- 
ern hemisphere. We expect that it is very unlikely that 
the spot and observer could be in opposite hemispheres, 
so we restrict our models to ones where 9 and i are both 
less than 90°. The parameters a and b are not freely 
varied: the anisotropy is restricted to < a < 1 and 
the parameter b is restricted to be within 1 sigma of the 
value given by the spectral model. 

In this paper we also perform simultaneous fits to the 
data from two time periods. In these simultaneous fits, 
the parameters M, R, and i are assumed to have the same 
value in both time periods. All other parameters are 
allowed to change their values, subject to the restrictions 
on 9, a and b mentioned in the previous paragraph. This 
brings the number of free parameters to 14. 

3.1. Effect of Spot Size, Shape and Number 

In this paper we restrict the models to only one visible 
spot. Prior to doing this we carried out fits to the data 
from all seven time periods with a 2-spot model. The 
best fit 2-spot models in all cases had zero amplitude 
for the second spot, indicating that no second spot is re- 
quired by the observed pulse sh apes. This is consistent 
with the model introduced by llbragimov fc Poutanenl 
(2009) where the accretion disk hides the antipodal spot 
on the other side of the star during the early stages 
of the outburst (stag es 3 and 4 in the terminology of 
lHartman et al.l (I2008D; stages P a nd SD in the termi- 
nology of llbragimov fc Poutanenl (|2009T > V In the later 
periods of the burst (sta ges 5 and 6) the disk recedes 
away from the star in the llbragimov fc Poutanenl |2009) 
model, exposing the antipodal spot and creating a more 
complicated pulse shape. The low flux and signal-to- 
noise in the pulse shapes corresponding to stages 5 and 
6 of the 2002 and 2005 outbursts mean that it is not 
feasible to use these p ulse shapes in the present analysis. 

In an earlier paper (jLeahv et al.ll2009l ). the light curve 
of a different neutron star, XTE 1814, was analyzed. In 
that case a one-spot model led to a residual (or bump) in 
the pulse shape that was seen in all wavelengths. In that 
case, the residual provided the motivation for a two-spot 
model for the data which did indeed fit the data signif- 
icantly better. In the case of SAX J1808, the one-spot 
models (shown in Section do not lead to a significant 
excess feature at all wavelengths that would motivate a 
t wo-spot model. 

iPoutanen fc Gicrlinski ([2003D have shown using ana- 
lytical approximations that the theoretical pulse shape 
resulting from a hot spot is independent of the spot 



size, as long as the spot radius is small compared to the 
size of the neutron star radius. Pre v ious computations 
of pu lse shapes (e.g. iLeahv et al.l (|2008l) iLeahv et all 
(2009)) have confirmed that the assumed spot size does 
not strongly affect the resulting pulse shape. For this 
reason, we use a simplified computational procedure that 
treats the spot size as infinitesimally small. This speeds 
up the computations considerably compared to the use 
of spots that cover a larger area. Since the distance to 
the star and the overall intensity normalization are kept 
arbitrary, the area of the spot is not needed in our mod- 
els. 

The best-fit models using one infinitesimal spot have 
unphysically low masses. For example in row 1 of Ta- 
ble [1] we show the best-fit model for the 2002B3 epoch 
with 2M/R fixed at 0.5. (Results for other values of the 
gravitational radius are similar.) Since there are 2 en- 
ergy bands with 32 time-bins each and a total of 8 free 
parameters, there are 56 degrees of freedom. The best fit 
model has a chi-squared value of 69.7 for 56 degrees of 
freedom, and a mass of only 0.3 Mq which is unreason- 
ably low for neutron stars, or for quark stars. Increasing 
the radius of the spot so that it has a radius of O.li? or 
0.2i? does not change the mass or \ 2 significantly, as can 
be seen in rows 2 and 3 of Table [T] 

We also tested the hypothesis that a non-circular spot 
could affect the results. We tested a long, thin uniform 
intensity spot that had a length of OAR azimuthally. The 
result for this extreme shaped spot is shown in row 4 of 
TableQ]where it can be seen that there is again no signifi- 
cant change from the infinitesimal spot. A similar test of 
a long, thin spot with only latitudinal extent gives simi- 
lar results. An alterate type of non-circular spot can be 
approximated by a two spot model where both spots are 
artifically forced to have the same intensity, but the rela- 
tive angular locations are allowed to be free. This would 
allow for two spots to be close together, giving an effec- 
tive irregular shape. The result for this two-spot model 
is shown in row 5 of Table [1] This model marginally im- 
proves the x 2 j but the best-fit mass is still only O.354M . 
We conclude that changes in the spot size and shape do 
not help the models to converge to physical values for 
neutron stars. 

3.2. Scattered- Light Model 

The hot-spot model described in section [3] was used to 
fit two different data sets simultaneously. The resulting 
fits were found to be very poor. For example, data from 
1998B4 and 2002B3 were fit using this model. The case 
for 2M/R — 0.5 is shown in the first row of Table [2] 
(labelled "No Scatter"). The resulting poor \ 2 of 200 for 
113 degrees of freedom is typical for fits done for other 
values of 2M /R and for other pairs of data. 

The poor fits that result from the application of the 
hot-spot model motivates us to introduce an extra com- 
ponent to the model that includes light that travels from 
the hot-spot on the star, scatters off of material near the 
star, and then arrives at the observer. The scattering 
material could be either an accretion column or an ac- 
cretion disk. If the light is scattered by the disk, then 
there will be a bright pattern on the disk that will appear 
to rotate around t he star at the same f r equen cy as the 
star, as studied by Sa zonov fc Sunvaevl (|2001l ). In this 
section we will outline a simple model for scattering off 
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of an optically thin cloud of electrons, such as the ionized 
surface layer of an accretion disk. 

In this simple model, we do not include the effects of 
gravitational light-bending or doppler boosting for the 
scattered light in order to include only the most impor- 
tant features. The reason why we can neglect these nor- 
mally important relativistic effects is that the amplitude 
of the scattered light is a very small fraction of the ampli- 
tude of the light that travels directly from the spot to the 
observer. We have tested that the inclusion of Doppler 
boosting to the scattered light model does not make any 
significant difference to the outcome of the fits. 

Consider a spot at co-latitude 9, azimuthal angle <f> on 
the surface of a star of radius R. The unit vector pointing 
from the center of the star to the spot is 



(cos 4>x + sin <py) + • 



(2) 



Photons emitted by the spot are scattered by material 
at co-latitude $ and azimuthal angle ip at a distance r. 
The unit vector pointing to the location of the scatterer 
is 

r = sin$(cos<^x + simpy) + cosi9z. (3) 

When light bending is neglected, the photon travels on 
a straight line defined by the vector £ = rf — Rh which 
has magnitude 



I = r (sin 2 a + (cos a — R/r) 



2\l/2 



(4) 



where a is the angle between the vectors f and n defined 

by 

cos rr = cos#cos-# + sin#sin$cos(0 — p). (5) 

The vector pointing towards the observer is k = sin ix + 
cosiz, so the photon must be scattered through an angle 
p in order to be detected, where p is defined by 



cosp = 



-l-k 



1 



[— i?(cos 9 cos i + sin 9 sin i cos ( 



+r (cos -d cos i + sin i) sin i cos p)] 



(6) 



The differential cross-section for Thompson scattering 
is g a (1 + cos 2 p) where p is the angle that the pho- 
ton is scattered through. If the flux of light originating 
from the hot spot at angle <fi that impacts the scatterer 
at p is F(<fi), the flux of scattered light reaching the ob- 
server is then F sc oc F((j>)(l + cos 2 p) where the constant 
of proportionality depends on the optical depth of the 
scattering material. 

In order to model scattering of light by an optically 
thin cloud of electrons close to the accretion disk, the 
co-latitude of the scattering material is set to # = tt/2. 
Photons emitted by the spot at angle <\> = 2irt/P + <f>o 
will impact the disk over a range of azimuthal angles ip 
and radii r. The pattern appearing on the disk will be 
centered around an angle tp equal to </>. However, the 
photons emitted by the spot and then scattered by the 
disk reach the observer at a later time than the photons 
that travel directly to the observer without scattering. A 
simple way to treat the time lag is to define an effective 
phase lag A<f> so that tp = <f> + A<fi, where Acj> includes 
the delays due to the light propagation time from the 
star to the disk as well as the light-crossing times for the 
photons to cross the disk. The lowest order expression 



for the phase lag (ignoring relativistic corrections) is 

a , „ ,r(cos 2 9 + (sm9 - R/r) 2 ) 1 / 2 r . . 

A6 = 2ir(-± i — '—L2 sin i cos 0) 

cP cP 

where P is the neutron star's spin period. The time lags 
for the scattered light should be much less than the star's 
spin period. Since cP = 750 km for SAX J1808.4-3658 , 
physical solutions should have a scattering radius r « 
cP. In this model, we are approximating the extended 
pattern on the disk as an effective point source. We 
tested more complicated models with a smeared out pat- 
tern on disk and found that the changes from the point 
source model were insignificant. 

The minimum radius of impact on the disk for a photon 
emitted at co-latitude 9 is r m i n — R/ sin9. Similarly, 
parts of the disk at radii r > r shadow can be observed 
at all values of tp, where the shadow radius is r s hadow — 
R/ cosi. As a result, the pattern on the disk is always 
visible to the observer if 9 < tt/2 — i. If t he spot is at 
a high latitude (ie., 9 <C 1) as suggested bv lLamb et al.l 
(2009), then the illuminated part of the disk will be at 
r» JJ and the illuminated part of the disk will not be 
eclipsed by the star. 

As the spot moves around the star, the resulting pat- 
tern on the disk appears to move around the star at the 
same rate. In the frame rotating with the star, the pat- 
tern will not appear to move, so the flux F(cf>) impacting 
the disk will be a constant. As a result, the scattered 
flux observed will be proportional to just the differential 
cross-section. In addition to the 14 parameters required 
for the hot spot model introduced in the previous sec- 
tion, the scattered light model requires that we add two 
more free parameters: I sc and r. The scattered light is 
added to the model by adding the term 



Fsc — -^sc(l 



COS' 



P) 



(8) 



to each energy band. Since we are implementing a simple 
model for scattering from a disk (■& = tt/2), the scattered 
flux is 



F sc = I sc \l+{j) [sinicos(0 + A(f>) 
R 

(cos 9 cos i + sin 9 sin i cos <fi) 

r 



0) 



where A<f> depends on the free parameters r, i, 9 through 
equations (0| and ([7]). The normalization for the hot spot 
model is defined so that the amplitudes I\ and li are of 
order unity, so a physical scattered light model should 
have an amplitude I sc << 1- 

In row 2 of table [2] the results of a joint fit to the 
1998B4 and 2002B3 data sets using the hot spot model 
along with the simple scattering model given by equation 
© are shown (labeled "model 1"). Three new parame- 
ters were introduced, the flux amplitudes I s \, I S 2 for the 
scattered light in 1998B4 and 2002B3 and r, the effec- 
tive radius at which the scattering mainly occurs. The 
addition of these 3 parameters to the model bring the % 2 
value down to 120, which is an acceptable fit for 110 de- 
grees of freedom. The resulting best-fit stellar model has 
a mass of 1.44M Q . The best- fit spot latitudes are close to 
the North pole in both time periods, while the inclination 
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angle is large, placing the observer close to the equatorial 
plane of the star. The best-fit scattering radius of If 6 km 
satisfies both inequalities r > Rj sin #1.2 and r > Rj cosi 
so that this location will be illuminated by the spot and 
will always be visible to the observer. The relative ampli- 
tude of the scattered radiation in f 998B4 is 0.2% of the 
flux coming directly from the spot. Similarly, the relative 
amplitude for the 2002B3 data is 0.8%. Both amplitudes 
are very small fractions of the direct flux, as would be 
expected for scattered light. The minimum in x 2 is fairly 
broad in the parameters corresponding the mass of the 
star and inclination, with the one-sigma allowed region 
including a range of masses between 1 .2 - 1 .6 M@. Only 
one model (for 2M/R = 0.5) is shown in tabled but the 
results are similar for other values of 2M /R. 

The best-fit solutions using the scattering model de- 
scribed above all have the property that the location of 
the scatterer is large compared to the radius of the star, 
r S> R. This motivates the approximation of equation 
([9]) to lowest order in R/r, 



F sc = I sc (f + sin 2 i cos 2 (0 + A</>)) 



(10) 



Furthermore, the inclination angles are fairly close to 90 
degrees, so sini ~ 1 . Using this approximation leads to 



F,, 



Isc (I 



\4> + &<t>)) 



(11) 



If the best- fit value of r/R 3> I and sini ~ 1, equation 
(fTT|) should be a good approximation to equation ([9]) . 

As an alternative to the scattered light model given 
by equation (j9), we now test the simpler model given 
by equation (|II[) . with A0 as a free parameter, instead 
of r. The best-fit stellar model that results when equa- 
tion (fTTj) is added to the basic hot-spot model and fit 
to the I998B4 and 2002B3 data sets is shown in row 3 
of table [2] (labelled "Model 2"). The best-fit star (with 
2M/R = 0.5) using Model 2 has x 2 = 120 for 1 1 degrees 
of freedom, so it is as good a fit as Model I. The best-fit 
values of the various parameters are different from the 
values for Model 1, however, the differences are smaller 
than the one-sigma differences allowed by either model. 
For instance the mass using Model 2 is 1.57 M Q , which 
lies in the one-sigma region allowed by Model I. We 
conclude that the two models produce results that are 
statistically the same. For this reason, we use the sim- 
pler model using equation ([TT]) to model the scattered 
light for the rest of this study. 

Since the addition of 3 parameters (J s i,2 and A<j)) 
brings about a significant reduction in x 2 ; and the masses 
are raised to realistic values, we include this scattered- 
light model in subsequent fitting of the pulse shapes. 
In summary, this component is required in order to ob- 
tain acceptable fits and to obtain physically reasonable 
masses. 

A small constant (or "DC") flux component was re- 
quired in a similar analysis of XTE J1814 (see section 
3.2 of iLeahv et all (|200l ). In the case of XTE JI814, 2 
free parameters (one for each energy band) were included 
to improve the quality of the fits. The amplitudes of the 
DC components were less than 2% of the pulsed compo- 
nent, similar to what is seen in the present analysis of 
SAX Jf 808. 

4. RESULTS 



A previous pulse-shape analysis ([Leahy et al.ll2008i ) of 
the 1998 outburst led to a very low best-fit mass and ra- 
dius for the neutron star. Although the 1998B4 data are 
taken o ver a slightly differe nt time period than the data 
used bv ILeahv et al.1 (|2008h . a reanalysis of the 1998B4 
data does not change the best fit mass and radius by very 
much. For instance, Table 5 of lLeahv et al.l ([2008) shows 
the best fit model has a neutron star mass of 0.96M Q and 
X 2 /dof = 36.3/24. When we fit the 1998B4 data includ- 
ing the scattered-light model, we find a best fit mass of 
0.98A/© and x 2 /dof = 63.3/54, which is not significantly 
different. 

An alternative strategy for fitting the data, which 
should reduce the effects of any near-degeneracies in 
model parameters, is to do a joint fit between the 1998B4 
data and the other data segments. As noted in section [2j 
the 1998B4 data has a much higher modulation than the 
other data segments. In addition, the error bars for the 
1998B4 data are relatively small and the data has fairly 
low scatter compared to most of the other data segments. 
In the x 2 comparison between the various data streams 
described in section[2]we note that the largest differences 
in pulse shape can be found when comparing the 1998B4 
data to any of the other data streams. For these reasons 
it makes sense to attempt to jointly fit the 1998B4 data 
with the other streams. 

For all joint fits we consider two data streams at a 
time, 1998B4 and one other data set. Due to the large 
number of parameters needed to model the data, it is not 
computationally feasible for us to jointly fit more than 2 
data streams at a time. In a joint fit, the mass, radius 
and inclination angle of the neutron star are assumed to 
be identical in all years. The angular location of the spot 
6, anisotropy parameter a, blackbody to Compton ratio 
b, scattered-light parameters, normalizations and phases 
are all allowed to vary with time. This leads to a total 
of 18 parameters in the joint fits. Since there are 4 light 
curves (2 energy bands for each time period) with 32 
points each, the total degrees of freedom are 110. 

4.1. Joint Fits to 1998B4 and 2002B3 Data 

The joint analysis of the 2002B3 data with the 1998B4 
data results in best-fit neutron star models with larger 
masses compared to the best-fit models that only use the 
1998B4 data. Table[3]shows the values of the best-fit pa- 
rameters for a range of neutron star compactness values 
(2M/R). Each row of Table[3]is generated by keeping the 
compactness ratio 2M/R (at the location of the spot) at 
a fixed value and varying all other parameters until the 
lowest value of x 2 is found. For each value of 2M/R the 
following variables are displayed: the mass of the neu- 
tron star M, the equatorial radius of the neutron star i?, 
the angular location of the spot 9\ (as measured from the 
North pole) in 1998, the angular location of the spot 82 in 
2002, the inclination angle i, the anisotropy parameters 
for 1998 (ai) and 2002 (02). We assume the orbital and 
spin axes are aligned, so i is also the binary's inclination 
angle. In addition there arc 4 normalizations, 2 phases, 2 
blackbody to Compton ratios, and 4 scattered-light pa- 
rameters which are not displayed in the table, since they 
are of no interest here. A x 2 penalty was used to keep the 
blackbody to Compton ratios b± and 62 within 1 sigma 
of the spectral models. 

The best-fit models in Table [3] should be compared 
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Figure 7. Mass and Radius confidence contours for joint fits to 
1998B4 and 2002B3 data using the constraint i < 90°. Mass- 
radius curves (solid curves) for stars spinning at a frequency of 
401 Hz are shown for a variety of E OS. Equatio ns of state dis- 
playe d are: Arnctt & Bowers ( 1977) mod els A l lPandharipandei 
1 19711) and L jPandharipandc ct al. 2006); mixed phase quark- 
hadron stars A BPR1- 3 (Alford ct al. 200% APR lAkmal et al.l 
119911); BBB1 IBaldo et all 1199^) ; hyperon stars H4 and H7 
(Lackcv. Navvar, & Owen 2006); and quark stars Q160 and Q180 
(where the number co rresponds to value of B 1 / 4 in MeV where 
B is the bag constant Glcndenning (2000)). Confidence contours 
for 2<t (dashed curve) and 3cr (dot-dashed curve) are shown. The 
boundaries shown as dotted lines correspond to the largest (0.6) 
and smallest (0.2) values of 2M/R used in the computations. 

with Table 5 of iLeahv et al.l (|2008|) where only data 
for 1998 was used. Consider the best-fit models with 
2M/R = 0.4 in both tables. When only 1998 data 
was used, the best-fit model had small mass and radius 
(Al = 0.90M Q and R = 6.7 km) and both the spot's an- 
gular location and the inclination angle were close to 30 
degrees. The addition of data from 2002B3 allows for a 
larger mass and radius (M = 1.43M Q and R = 10.8 km) 
for the same value of compactness. The larger radius is 
compensated for by moving the spot closer to the spin 
axis with B\ close to 10 degrees in 1998 and 5 degrees in 
2002B3. 

The best-fit values of the scattered- light amplitudes A\ 
are 2 x 10~ 3 for 1998B4 and 8 x 10~ 3 for 2002B3, where 
the flux is normalized to 1. 

The best-fit model with 2M/R = 0.6 in Table [3] is 
shown in Figures [1] and [5] with dashed curves for the 
low-energy bands and solid curves for the high-energy 
bands. All four light curves are fit simultaneously, but 
for clarity, the data and best-fit models are plotted in 
two separate figures. The relative phases between the 
high and low energy-bands for each data set are included 
in the models. 

Figure [7] shows 2 and 3cr confidence contours on the 
mass-radius plane for joint fits with 1998B4 and 2002B3. 
This figure is generated by specifying a mass and then 
varying all parameters until the lowest value of % 2 for 
the mass is found. Two constraints were added to the 
variation of parameters. First, the ratio 2M/R was only 
sampled in the range 0.2 < 2M/R < 0.6 in order to stay 
within the bounds of reasonable equations of state. Sec- 
ond, an unconstrained variation of the inclination angle 
leads to acceptable fits for solutions with i > 90°, which 
is very unlikely. Hence, we constrained the inclination 
angle to values with i < 90°. The resulting confidence 
contours in Figure [7] should be compared with the cor- 
responding confidence contours for the 1998 data shown 
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Figure 8. Inclination ranges for the joint fit to 1998B4 and 
2002B3 data sets. Points labeled with a triangle correspond to 
the best-fit models shown in Table [3] For each one of these points, 
the 3<r range of inclination angles (in degrees) for the given mass 
and radius is shown. The dot-dashed curve shows the same 3o" 
limit on the allowed values of mass and radius shown in Figure [7] 
Selected points along the 3cr contour are shown with circles. The 
inclination angle corresponding for each of these points is shown. 

with bold curves in Figure 3 of ILeahv et all (l200l (la- 
belled "oblate and time delays"). The 3<r confidence re- 
gion only allows very small stars with R < 7 km if only 
the 1998 epoch is used. The additional data from 2002B3 
and allowance for scattered-light used in the present anal- 
ysis allows stars with larger radii, extending to 13 km at 
the 3(7 confidence level. The largest mass stars allowed 
(at 3<r) have M = 1.7M Q . At 3c, the data allow for a 
wide variety of modern equations of state. 

4.2. Inclination Constraints 

The best-fit models displayed in Table [3] have inclina- 
tions ranging from 70 - 90 degrees. However, for each 
best-fit set of values for mass and radius shown in Ta- 
ble El there is a range of inclination angles allowed by the 
data. The range of inclination angles allowed at the 3 a 
level is found by keeping the mass and radius fixed at the 
values shown in the table, and allowing all other parame- 
ters to vary. In FigurelH the locations on the mass-radius 
plane for the five models shown in Table [3] are marked 
with triangles. The allowed 3er range of inclination an- 
gles (in degrees) for each of these points is shown on 
Figure [5] For example, the model with M — 1.43M 
and R — 10.8 km is allowed inclination angles in the 
range of 73° < i < 81°, although the lowest \ 2 results 
for i = 75.5°. 

On Figure [8] the same 3 a contours on the mass-radius 
plane that were shown in Figure [7] are redisplayed. For 
selected points on the contours (shown with circles), the 
one value of inclination that is allowed is shown. For 
example, Figure [5] shows that the data allows at the 3 
a confidence level a star with M = 1.4AM®, R = 12.5 
km and i = 79°. In general, the lowest mass stars re- 
quire smaller inclination angles (near 40 degrees) while 
the highest mass stars require inclinations that are as 
high as is allowed. 

There are a number of observations that suggest that 
the inclination angle may be smaller than 90 degrees. 
The strongest evidence is the absence of eclipses, which 
cons trains i < 82° assuming a Roch e lobe filling compan - 
ion (jChakrabartv & Morgan! I1998D . I Wang et al.l (2001) 
observed V4580 Sgr, the optical counterpart of SAX 
J1808, during the 1998 outburst. They modelled the op- 
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Figure 9. Mass and Radius confidence contours for joint fits to 
1998B4 and 2002B3 data using the constraint i < 70°. Mass- 
radius curves (solid curves) for stars spinning at a frequency of 401 
Hz are displayed with solid curves. The radius corresponds to the 
equatorial radius of the star. EOS labels are the same as in Figure 

Confidence contours for 2 a (dashed curve) and 3 a (dot-dashed 
curve) are shown. 

tical and IR data using an X-ray-heated accretion disk 
and found limits on the inclination of 42° < i < 88° 
if the distance is d = 2 kpc and 20° < i < 65° if 
d = 3 kpc. Further model ling of the optical data ob- 
served during quiessence bv lWang et all (12009) suggests 
that i < 70° using the distan ce of d = 3.5 kpc deri ved by 
iGallowav k Cummingj pOOl ). iDelove et all ((2008) mod- 
elled the optical data from quiessence using a model that 
includes emission from bot h the disk and the c ompanion. 
Using a distance of 3.5 kpc. lDelove et all (|2008T ) constrain 
the inclination to 32° < i < 74° if 10% uncertainty in the 
distance is assumed. 

A feature interpreted as an iron line relativistically 
broadened through the orbital motion of the disk was 
de tected in X-ray observatio ns during the 2 008 outburst 
bv lCackett et all (|2009D and lPapitto et al.l (|2009fl . Both 
groups modelled the emission feature and constrained the 
inclination angle. iCackett et al.l (|2009T ) found i = 55+® 
degrees, while IPapitto et all ([2009D provide a less precise 
constraint of approximately i > 60°. A different type 
of constraint on the inclination arises from the analy- 
sis of the 2002 outburst data by llbragimov fc Poutanenl 
( 2009). In their analysis, llbragimov fc Poutanenl (|2009h 
assume that the X-ray emission is produced by two an- 
tipodal spots, one of which is sometimes hidden by the 
accretion disk, due to movement of the inner edge of the 
accretion disk. In this model, they find that the data can 
be explained by inclinations in the range i < 73°. 

Although all of the constraints on inclination have 
some model dependence, there is a general agreement 
that inclinations smaller than about 70° are suggested 
by the different data. With this in mind, we refit the 
data for 1998B4 and 2002B3 subject to the constraint 

1 < 70°. The best fit models for fixed values of 2M/R 
are shown in Table 2] Comparing Table 0] with Table 
[3] it can be seen the the most significant change is that 
lower mass stars are allowed when the inclination angle 
is constrained to i < 70°. The values of x 2 increase, but 
the magnitude of the increase is not significant. 

Figure [5] shows the 2 and 3-er confidence contour on 
the mass-radius plane using the constraint i < 70° for 
joint fits to the 1998B4 and 2002B3 data. The 3 a confi- 
dence region shrinks when the constraint on inclination 



is added, however stars as massive as 1.5M© are still 
allowed. 

4.3. Location of the Hot Spot 

In all of the hot spot models for the 1998B4 and 2002B3 
data sets, the colatitude of the spot 9 is larger (further 
from the pole) in 1998 compared to the colatitude in 
2002. This trend is seen in all models fitting the other 
data sets as well. Since a hot spot located close to the 
spin axis produces a less-modulated pulse shape, this is a 
natural outcome of the fact that the 1998 data is highly 
modulated (~ 10%) and the data for the later outbursts 
have smaller modulation (~ 5%). 

In the best-fit models for the 1998B4 and 2002B3 data 
shown in Table [3] the shift in colatitude A9 = 9\ — 62 
ranges from 5° to 8°. When all models that fit the data 
within 3 a confidence are considered, the colatitude in 
1998 lies in the range 9° < 9 X < 32°, the colatitude in 
2002B3 lies in the range 4° < 9 2 < 17°. The 3 a range 
of colatitude difference is in the range 6° < Ad < 15°. 
This wandering of the h ot spot locati o n is co nsistent with 
the model proposed by lLamb et al.l (|2009f) where a hot 
spot close to the spin axis moves due to movement of the 
star's magnetic field. 

4.4. Joint Fits to 1998B4 and 2002B4 Data 

The 2002B4 and 1998B4 data have similar pulse shapes 
(although a statistical test shows that they are different 
at high significance - chance probability < 10~ 20 ). Of all 
the data sets, these two epochs have the smallest error- 
bars which are roughly equal in magnitude. The most 
important difference is in the modulation of the pulse 
shapes. 

The best-fit models for joint fits to the 1998B4 and 
2002B4 data sets are shown in Table [SJ In this ta- 
ble, the subscript "2" corresponds to values of in the 
2002B4 epoch, while the subscript "1" corresponds to 
the 1998B4 epoch. These best-fit models should be com- 
pared with the best-fit models shown in Table [3l The 
range of masses and radii found when the 2002B4 data 
are fit is similar to those found when data from 2002B3 
are used. The values of mass and radius shown in Table[5] 
are consistent with the 3cr limits arising from the joint fit 
between the 1998B4 and 2002B3 data, with the excep- 
tion of the model with a mass of 0.66Af Q . The best- fit 
values of the 1998 angular location of the spot 9\ and the 
inclination angle i differ in the best-fit models displayed 
in Tables [3] and [5j This is not surprising, since as we 
have shown earlier, given values of mass and radius, a 
wide range of angles 9 and i will give acceptable fits to 
the data. 

The best-fit model's light curve for the 2002B4 data is 
shown in Figure [3l (The light curves for 1998B4 gener- 
ated by the best-fit model are not distinguishable from 
the curves shown in Figure [T] and are not plotted.) 

4.5. 2005 Outburst Data 

The 2005 outburst was only about 70% as luminous 
as the 2002 outburst (|Hartman et al.ll2008| ) . In addition, 
the RXTE satellite has lost sensitivity since 2002. As a 
result, the 2005 data is much noisier and has larger error 
bars than the 1998 and 2002 data as can be seen in Fig- 
ures |4]|6l One-spot models were used to simultaneously 
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fit the 2005 data with the 1998 data. The best-fit model 
parameters are shown in Tables [6][8] and the best-fit light 
curves are shown in Figures [4][6j 

The best-fit models that make use of the 2005 outburst 
data shown in Tables [6] -[8] are very similar to the best-fit 
models computed using the 2002B4 data and displayed 
in Table 03 The models with 2M/R = 0.4,0.5 and 0.6 
(and with M > 1.0M©) are consistent with the 3<r limits 
arising from the joint 1998B4-2002B3 fits. 

The 2008 outburst was almost as bright as the 2005 
outburst, however the degradation of the RXTE PCA 
over the intervening 3 years has t he result that the puls e 
shapes from the 2008 outburst (|Hartman et al.1 120091) . 
while similar to the 2005 pulse shapes, are even noisier 
than the 2005 pulse profiles. For this reason, we have 
not attempted to fit the 2008 pulse profiles. 

5. DISCUSSION 

In this paper we analyze multiple pulse shapes result- 
ing from the different epochs of the 1998, 2002 and 2005 
outbursts of SAX J1808. We jointly fit data from 1998 
and several other epochs in order to find consistent mod- 
els for the neutron star's mass, radius and inclination. 

In each joint fit to the data from two different epochs, 
the mass, equatorial radius and the inclination angle 
are kept fixed between epochs, while the spot loca- 
tion, anisotropy and blackbody-to-compton ratios are al- 
lowed to change from epoc h to epoch. We use a spec- 
tral model, motivated by [G ierlins ki et al.l <|2002f ) and 
iPoutanen fc Gierlifiskil (|2003[ ). that includes blackbody 
radiation and Comptonized radiation in the low-energy 
(3-4 keV) band and only Comptonized radiation in the 
high-energy (9-20 keV) band. We find that an extra 
scattered-light component is required to fit all of the 
pulse shapes. 

The inclusion of the scattered radiation model (Sec- 
tion [XU, together with a hot spot that wanders a little 
in its co-latitude, allows a consistent fit for all different 
epochs with a consistent mass. Increased radius or in- 
creased spot co-latitude increases the second harmonic 
content in the pulse shape; scattered radiation reduces 
it; and increased light bending (smaller star at a given 
mass) can increase or decrease second harmonic content 
depending on viewing geometry. Without including all 
these physical effects, the fits to some epochs require a 
small star whereas other epochs require a larger star. As 
we have shown, by including all effects, we can consis- 
tently fit different epochs with a consistent mass. 

The joint fit of data from the 1998B4 and 2002B3 
epochs using the constraint i < 90° leads to the 3a con- 
fidence constraint on the neutron star mass 0.81M Q < 
M < 1.72M and equatorial radius 5.0km < R < 
12.7km, as summarized in Figure [71 Inclinations as low 
as 41° are allowed. If the inclination is further con- 
strained to i < 70°, as suggested by many separate ob- 
servations, the 3cr limits are 0.83M Q < M < 1.5OM 
and 7km < R < 13km, as summarized in Figure [9j In- 
clinations as low as 35° are allowed in this case. 

The joint fits between 1998B4 and the other epochs 
(2002B4, 2005B1, 2005B2 and 2005B4) yield indepen- 
dent best-fit models summarized in Tables [SJSJ The best- 
fit models from these epochs with M > I.OMq have mass 
and radius values that lie within the 3er confidence region 
for the 1998B4-2002B3 joint fits. This indicates that the 



various joint fits between the different epochs give con- 
sistent results. 

The 3cr confidence region for the 1998B4-2002B3 joint 
fits allows a much wider range of masses and radii for 
the neutron star than a simil ar fit that only mad e use of 
data from the 1998 o utburst dLeahv et all [20081 ). In the 
previous analysis by iLeahv et al.l (|2008f ) a similar type 
of analysis using 2 narrow energy bands constrained the 
neutron star mass and radius to very small values with 
R < 7 km. The inclusion of a scattered-light component 
and use of joint fits with data from other epochs (2002 
and 2005 outbursts) results in stars with larger masses 
and radii than are given by fits to the 1998 data alone 
without the scattered-light component. 

In their analys is of the 2002 data, 

libra gimov fe Poutanenl (|2009f ) do not attempt to 
find a best-fit mass and radius for SAX J1808. For some 
of their models they make use of a fixed mass of 1.4 
Mq and a radius of either 10 or 12 km. However, the 
results of their modelling (such as the magnetic moment, 
visibility of an antipodal spot) do not strongly depend 
on the assumed values of mass and radius. A mass of 
1.4 Mq and a radius in the 10-12 km range is consistent 
with our final results. 

Observatio ns of SAX J1808 in quiescence by 

iHeinke e~ al. (2007. [20091) have set limits on the luminos- 
ity. The low luminosity during quiescence implies that 
rapid neutrino cooling takes place in the inner core, such 
as the direct URCA process. lYakovlev fe Pethickl (|2004[ ) 
show that direct URCA can take place if the core density 
is larger than approximately 10 15 g/cm 3 . For a stan - 
dard hadronic EOS such as APR i|Akmal et aT1ll998[) , 
stars with masses larger than about 1.5M® (and spin- 
ning at 400 Hz) have central densities larger than this 
critical density. Softer EOS will have lower mass stars 
for the same central density. This suggests that many of 
the models allowed by the 3er confidence region shown in 
Figure [7] could have cores dense enough to allow cooling 
through the direct URCA process. 

A pulse-shape analysis o f the accreting ms -period X- 
ray pulsar XTE 1814-338 (|Leahv et all [2009) using the 
similar methods as used in this paper implied that the 
mass and radius of XTE 1814-338 must be quite large. 
Although the 3cr regions for XTE 1814-338 and SAX 
J1808 do not overlap, the results are not necessarily in- 
consistent. It is only necessary for an EOS mass-radius 
curve to pass through the regions for both stars. For 
example, an EOS mass-radius curve that is slightly to 
the right of APR in Figure [7] would still be allowed by 
the SAX J1808 data, and would be stiff enough to be 
allowed by the XTE 1814 data, and also allowed by the 
const raints on the slow X-ray pulsar Hercules X-l ([Leahvl 
12001 . 
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Table 1 

Effect of Spot Size and Shape for stellar models with 2M/R = 0.50 fit with data from 2002B3. 



Model 


M 


X 2 /dof 


i 


a 


02 


64 




Mq 


deg. 


deg. 








infinitesimal spot 


0.304 


69.7/56 27.0 


36.0 


0.50 






circle radius 0.1 R 


0.305 


69.7/56 27.0 


36.2 


0.51 






circle radius 0.2R 


0.307 


69.7/56 27.2 


36.6 


0.51 






azimuthal line OAR 


0.310 


69.2/56 26.2 


35.8 


0.51 






two spots (equal amp/ 


1 0.354 


62.1/54 12.0 


26.0 


0.54 


54.4 


60.1 






Table 2 










Best Fit Parameters for Joint Fits using Data from 1998B4 and 2002B3 with and without the scattered light component for 






2M/R = 0.5. 










M 


R 6»i 




0,2 Isl 


Is2 


r 


X 2 /dof 


Mq 


km deg. 


deg. deg. 






km 




No Scatter 1.02 


6.0 22.6 


10.6 49.5 0.57 


0.49 







200/113 


ivioaei i 1.44 


8.5 12.0 


5.6 80.0 0.80 


0.68 0.002 


0.008 


110 


i on /i in 
izU / iiu 


Model 2 1.57 


9.4 10.5 


4.8 88.9 0.92 


0.77 0.003 


0.007 




120/110 


Note. — Subscript "1" 


corresponds to 1998B4 values and "2" to 2002B3. 














Table 3 












Best Fit Parameters for Joint Fits using Data from 1998B4 and 2002B3. 






2M/R M 


R 


01 02 


i 


ffli 


0,2 


X 2 /dof 


M @ 


km 


deg. deg. 


deg. 








0.60 1.31 


6.5 


14.3 6.2 


89.6 


0.77 


0.62 


116/110 


0.50 1.57 


9.4 


10.5 4.8 


88.9 


0.92 


0.77 


120/110 


0.40 1.43 


10.8 


10.4 4.9 


75.5 


0.85 


0.73 


123/110 


0.30 1.19 


12.1 


10.0 4.7 


70.8 


0.89 


0.75 


124/110 


0.20 0.80 


12.5 


10.2 5.0 


70.3 


1.0 


0.88 


123/110 


Note. — Subscript "1" 


corresponds to 1998B4 values and "2" to 2002B3. 














Table 4 










Best Fit Parameters for Joint Fits using Data from 1998B4 and 2002B3 with i 


< 70°. 




2M/R M 


R 


01 02 


i 


ai 


a 2 


XVdof 


Mq 


km 


deg. deg. 


deg. 








0.60 1.21 


6.0 


18.4 8.3 


70.0 


0.64 


0.52 


118/110 


0.50 1.27 


7.6 


16.0 7.6 


63.3 


0.65 


0.56 


122/110 


0.40 1.30 


9.7 


12.7 6.1 


64.8 


0.72 


0.62 


124/110 


0.30 1.09 


11.0 


11.9 5.8 


61.9 


0.75 


0.65 


125/110 


0.20 0.80 


12.5 


10.2 5.0 


70.0 


0.99 


0.87 


123/110 



Note. - 



Subscript "1" corresponds to 1998B4 values and "2" to 2002B3. 
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Table 5 

Best Fit Parameters for Joint Fits using Data from 1998B4 and 2002B4. 



2M/R 


M 


R 


0i 


02 


i 


ai 


a-2 


X 2 /dof 




Mq 


km 


deg. 


deg. 


deg. 








0.60 


1.27 


6.3 


22.3 


18.6 


53.1 


0.57 


0.65 


105/110 


0.50 


1.27 


7.5 


21.0 


17.5 


46.9 


0.57 


0.64 


106/110 


0.40 


1.19 


8.9 


20.2 


16.6 


41.9 


0.57 


0.64 


111/110 


0.30 


0.96 


9.7 


19.5 


16.3 


39.4 


0.57 


0.64 


108/110 


0.20 


0.66 


10.1 


19.2 


16.4 


38.3 


0.58 


0.65 


109/110 


Note. - 


Subscript "1" 


corresponds to 1998B4 values and 


"2" to 2002B4. 


















Table 6 














Best Fit Parameters for Joint Fits using 


Data from 1998B4 and 2005B1. 






2M/R 


M 


R 


0i 


02 


i 


ffli 


a 2 


X 2 /dof 




Mq 


km 


deg. 


deg. 


deg. 








0.60 


1.36 


6.7 


19.6 


19.2 


57.2 


0.58 


0.66 


111/110 


0.50 


1.33 


7.9 


18.8 


18.6 


49.9 


0.58 


0.65 


114/110 


0.40 


1.20 


9.0 


18.3 


18.1 


44.8 


0.58 


0.65 


116/110 


0.30 


0.98 


9.9 


18.0 


17.9 


41.6 


0.58 


0.65 


117/110 


0.20 


0.67 


10.3 


17.9 


17.9 


40.1 


0.58 


0.66 


119/110 


Note. - 


Subscript "1" 


corresponds to 199SB4 values and 


"2" to 2005B1. 


















Table 7 














Best Fit Parameters for Joint Fits using 


Data from 1998B4 and 2005B2. 






2M/R 


M 


R 


0i 


02 


i 


ai 


a-2 


X 2 /dof 




Mq 


km 


deg. 


deg. 


deg. 








0.60 


1.20 


5.9 


23.1 


13.9 


53.2 


0.57 


0.64 


109/110 


0.50 


1.20 


7.1 


21.9 


13.1 


46.6 


0.57 


0.63 


109/110 


0.40 


1.09 


8.1 


21.1 


12.5 


42.1 


0.57 


0.62 


111/110 


0.30 


0.90 


9.1 


20.5 


12.3 


39.0 


0.57 


0.62 


113/110 


0.20 


0.64 


9.7 


20.1 


12.0 


37.5 


0.58 


0.63 


114/110 


Note. - 


Subscript "1" 


corresponds to 1998B4 values and 


"2" to 2005B2. 


















Table 8 














Best Fit Parameters for Joint Fits using 


Data from 1998B4 and 2005B4. 






2M/R 


M 


R 


0i 


02 


i 


ai 


a 2 


X 2 /dof 




Mq 


km 


deg. 


deg. 


deg. 








0.60 


1.19 


5.9 


23.0 


11.3 


54.4 


0.57 


0.64 


109/110 


0.50 


1.18 


7.0 


21.9 


10.7 


47.3 


0.57 


0.63 


111/110 


0.40 


1.07 


8.0 


21.1 


10.3 


42.5 


0.57 


0.63 


112/110 


0.30 


0.89 


8.9 


20.3 


9.9 


40.0 


0.58 


0.63 


114/110 


0.20 


0.63 


9.6 


19.9 


9.7 


38.3 


0.58 


0.64 


115/110 



Note. - 



Subscript "1" corresponds to 1998B4 values and "2" to 2005B4. 



